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Summary  (UL) 

^■Results  are  described  of  a  numerical  study  of  the  initial 
compressive  buckling  and  post-buckling  behaviour  of  GRP  panels 
reinforced  by  longitudinal  hat-section  stiffeners.  Particular 
reference  is  made  to  the  effects  of  interaction  between  local 
buckling  of  the  panel  laminate  and  overall,  column-like  buckling 
involving  bending  of  the  stiffeners.  Conclusions  are  reached 
regarding  the  influence  of  interactive  effects  on  collapse  strength 
and  recommendations  are  made  on  allowance  for  these  effects  in 
design. 
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INTERACTIVE  BUCKLING  EFFECTS  IN  STIFFENED  FRP  PANELS 


INTRODUCTION 

1.  Stiffened  GRP  panels,  employing  hat-section  stiffeners  aligned 
longitudinally  in  the  direction  of  dominant  load,  have  been  adopted 
in  the  deck  and  bottom-shell  structure  of  new  minehunters  currently 
under  development  for  the  Royal  Navy  (1)*.  A  theoretical  and 
experimental  investigation  (2)  has  identified  various  forms  of 
instability  which  may  occur  in  panels  of  this  type  under  longitudi¬ 
nal  compression  and  has  examined  the  relationship  between  initial 
buckling  modes  and  post-buckling  and  collapse  behaviour.  It 
appeared  that  significant  interactions  could  occur  in  some  cases 
between  local  and  overall  buckling  modes,  influencing  both  initial 
buckling  and  post-buckling  behaviour  and  hence  possibly  requiring 
special  treatment  in  design.  The  purposo  of  the  present  report  is 
to  describe  a  further  study  in  which  such  interactive  effects  were 
examined  by  numerical  analysis  of  some  panels  in  which  the  relation¬ 
ship  between  local  and  overall,  column-like  buckling  was  varied 
systematically  by  adjusting  the  panel  geometry. 

METHODS  OF  ANALYSIS 

2.  Methods  of  evaluating  the  initial  buckling  behaviour  of  stiff¬ 
ened  FRP  panels  have  been  described  previously  (2)  but  will  be 
summarised  below  for  completeness. 

3.  A  typical  GRP  deck  or  shell  panel,  reinforced  by  longitudinal 
stiffeners  of  trapezoidal,  top-hat  section,  is  shown  in  Figure  1. 
Strips  of  laminate  forming  the  panel  cross-section  may  be  related 
to  x  and  y-axes  parallel  to  their  longitudinal  and  transverse 
edges.  For  most  marine-type  GRP  laminates  based  on  chopped-strand , 
woven-roving  or  unidirectional  reinforcement  it  is  possible  to 
assume  that  each  strip  of  laminate  is  orthotropic  with  principal 
directions  of  elasticity  parallel  to  its  edges  and  that  plies  of 
reinforcement  are  stacked  symmetrically  about  mid-thickness  so  that 
no  direct  coupling  occurs  between  membrane  and  bending  deformations 
In-plane  deformations  of  each  strip  then  satisfy  the  orthotropic 
plane  stress-strain  relationship 


where  t  ,  E  ,  G  ,  y  and  y  are  young's  moduli,  shear  modulus  and 

M  J  *7  *  / 

Poisson  ratios  referring  to  the  principal  directions.  The  moment- 
curvature  relationship  for  the  strip  is 


*(  )  References  on  Page  11 
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are  flexural  rigidities  per  unit  width  in  x  and  y-directions. 

D  -  G  h3/12  is  the  twisting  rigidity  and  h  is  the  laminate 
thickness. 

I 

4.  In  the  presence  of  a  uniform  destabilising  stress  »x  as  shown 

in  Figure  1«  the  equation  governing  small  bending  displacements  w 
of  a  laminate  strip  is 


°x  ^ 
ax’ 


♦  2h 


ax  ay 


♦  D„ 


a4w 


N 


a2w 

*17 


13) 


where  H  *  2Dxy  +  «XDX  and  Nx  ■  «xh.  A  conservative  estimate  of 

local  buckling  stress  for  a  stiffened  panel  may  be  obtained  by 
applying  to  each  element  of  the  cross-section  the  formula 


N 


xcr 


/w 


(4) 


which  corresponds  to  a  solution  of  equation  3  for  a  long  ortho¬ 
tropic  strip  of  width  b  with  simply  supported  edgesi  local  buckling 
will  tend  to  occur  predominantly  in  the  most  rlender  element  with 
a  number  of  half-waves  over  the  length  (a)  of  the  panel  given 
approximately  by 


n 


15] 
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A  corresponding  upper-bound  estimate  of  initial  buckling  stress  is 
provided  by  the  formula  for  a  long  orthotropic  strip  with  clamped 
edges  (3) 


xcr 


^(2.4H  +  4.6 
b2 


16] 


in  which  case  the  preferred  number  of  buckling  half-waves  is 
approximately 


n 


17] 


5.  Initial  buckling  may  also  take  the  form  of  column-like  insta¬ 
bility  of  longitudinal  stiffeners  and  attached  shell  laminate 
between  heavy  transverse  frames  or  bulkheads.  Ignoring  rotational 
restraint  provided  by  the  transverse  members,  this  buckling  stress 
may  be  estimated  using  the  Euler  formula 


in  which  El  is  the  flexural  rigidity  of  a  stiffener  with  assumed 
effective  breadth  of  shell,  A  is  the  total  cross  section  of  the 
stiffener  with  attached  strip  of  shell,  a  is  the  spacing  of 
transverse  frames  or  bulkheads  and  GA#  is  the  shear  rigidity  in 

which  A(  may  be  taken  as  the  area  of  the  stiffener  webs. 

6.  More  accurate  (theoretically  exact)  evaluation  of  initial 
buckling  stresses  may  be  achieved  using  folded-plate  analysis,  in 
which  out-of-plane  displacements  of  each  strip  satisfy  equation  3 
while  in-plane  displacements  u  and  v  satisfy  the  simultaneous 
equations 
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If  conditions  of  simple  support  are  assumed  at  the  ends  of  each 

2  2 

strip  (v  ■  3u/3x  ■  0  and  w  ■  3  w/3x  *  0  at  x  *  o,  a)  solutions  of 
equations  3,  9  and  10  are  periodic  in  x  and  may  be  obtained  in  a 
way  which  accounts  rigorously  for  interactions  between  strips  of 
laminate  forming  a  stiffened  panel.  Details  of  this  analysis 
method ,  together  with  generalized  computer  programs  developed 
independently  by  several  investigators,  are  contained  in  References 
4  to  8. 

7.  Evaluation  of  nonlinear  and  post-buckling  behaviour,  including 
the  effects  of  imperfections,  may  be  carried  out  using  finite 
element  analysis.  Results  described  in  the  present  paper  were 
obtained  using  the  general  purpose  nonlinear  FE  Program  ASAS-NL, 
developed  by  Atkins  RtD  in  collaboration  with  ARE,  Dunfermline  (9). 
Eight-node  isoparametric  quadrilateral  shell  elements  were  employed 
with  constant  anisotropic  material  properties  as  defined  in 
Figure  2:  large  displacements  were  represented  using  an  updated 
Lagrangean  (moving  coordinate)  procedure  with  application  of  a 
modified  Newton-Raphson  iterative  equilibrium  correction.  The 
adequacy  of  FE  models  was  tested  by  obtaining  estimates  of  initial 
buckling  stresses  and  modes  by  linear  eigenvalue  analysis  for 
comparison  with  exact  folded-plate  results. 

SCOPE  OF  ANALYSIS 

8.  As  a  basis  for  evaluating  interactions  between  local  buckling 
and  overall  column-like  buckling,  analysis  has  been  carried  out  for 
a  series  of  longitudinally  stiffened  panels  having  the  same  cross- 
sectional  geometry  and  material  properties  as  the  experimental 
panel  L2  described  in  Reference  2.  By  varying  the  span  a  between 
transverse  frames  a  range  of  structures  was  obtained  in  which 
initial  column  buckling  stresses  varied  from  less  than,  to  substan¬ 
tially  more  than,  local  laminate  buckling  stresses.  Panel  dimen¬ 
sions  and  material  properties  are  indicated  in  Figure  2.  In 
addition  to  the  experimental  length  (a  »  3060  mm)  spans  of  1000, 
1500,  2000,  23S0,  2670  and  3500  mm  were  considered. 

9.  Initial  buckling  stresses  were  estimated  approximately  using 
equations  4  to  8,  with  the  breadth  b  taken  equal  to  the  unsupported 
span  between  stiffener  webs,  and  were  computed  accurately  using 
folded-plate  and  finite  element  analysis.  Stiffener  flanges,  which 
have  been  shown  to  influence  local  buckling  significantly  (2),  were 
included  in  folded-plate  and  finite  element  models.  Finite  element 
analysis  was  also  used  to  examine  the  nonlinear,  post-buckling 
behaviour  of  panels  including  the  influence  of  assumed  imperfec¬ 
tions.  As  shown  in  Figure  3,  the  finite  element  model  referred  to 
a  single  longitudinal  stiffener  with  attached  shell  laminate 
extending  over  two  half-spans  a/2.  Conditions  of  simple  support 
were  assumed  at  the  position  of  the  transverse  frame  EF,  with  a 
transverse  plane  of  symmetry  at  mid-span  position  CD  and  a  "moving" 
plane  of  symmetry  at  mid-span  position  ABt  in  this  way  account  was 
taken  of  "continuous  beam"  interactions  between  adjoining  spans. 
Longitudinal  planes  of  symmetry  AC  and  BD  were  assumed  to  occur 
mid-way  between  longitudinal  stiffeners. 
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10.  For  the  purpose  of  nonlinear  analysis,  initial  deformation  of 
each  panel  was  assumed  to  have  the  form  of  the  preferred  local 
buckling  mode  with  amplitude  wQp,  combined  with  an  overall  distor¬ 
tion  of  amplitude  wQg  having  the  form  of  the  interframe  column 

buckling  mode.  In  each  case  two  levels  of  imperfection  were 
considered: 

(i)  wQp/a  *  0.0005,  wQg/a  =  0.0013  (corresponding 

approximately  to  maximum  distortions  measured  in 
experimental  Panel  L2); 

(ii)  wop/a  =  0.00005,  wQg/a  =  0.00013. 

Initial  deformation  was  taken  to  be  antisymmetric  about  EF, 
positive  (towards  the  stiffener)  at  AB  and  negative  at  CD. 

Analysis  was  carried  out  by  incremental  application  of  uniform 
end-shortening  displacement  on  the  plane  AB  of  each  panel. 

DISCUSSION  OF  RESULTS 


Initial  Buckling  Behaviour 

11.  Initial  buckling  stresses  obtained  from  the  approximate 
formulae  and  by  folded-plate  and  FE  analysis  are  listed  in  Table  1. 
A  selection  of  corresponding  buckling  modes,  obtained  from  the 
folded-plate  solutions,  is  shown  in  Figure  4.  Initial  buckling 
stresses  determined  by  FE  eigenvalue  analysis  are  generally  within 
5%  of  exact  folded-plate  values,  lending  some  credibility  to  the  FE 
model  used  in  the  nonlinear  analysis. 

12.  As  observed  previously,  critical  stresses  corresponding  to 
local  buckling  of  the  shell  laminate  between  stiffeners  occur  in 
closely  spaced  pairs;  asymmetric  modes  (Type  B)  are  usually  (but 
not  always)  preceded  by  symmetric  (Type  A)  modes.  The  latter 

.  generally  involve  virtually  (but  not  exactly)  zero  displacement  of 
stiffeners  out  of  the  panel  plane.  Type  B  modes  involve  distortion 
of  stiffener  cross-sections  including  some  sideways  bending  of  the 
stiffener  tables.  As  might  be  expected,  local  buckling  stresses 
and  associated  modal  wavelengths  are  insensitive  to  the  span  a. 

13.  In  the  longest  panel  examined  (a  **  3500  mm)  overall  buckling 
occurs  with  virtually  no  local  deformation  at  a  critical  stress  6% 
less  than  that  given  by  the  Euler  formula  (Equation.  8).  As  panel 
length  is  reduced  the  "column”  mode  becomes  increasingly  influenced 
by  local  deformation  of  the  shell  laminate  and  the  associated 
critical  stress  falls  progressively  below  the  Euler  value.  In  the 
shortest  panel  examined  (a  *  1000  mm)  local  shell  deformation 
dominates  the  symmetric  n  •  1  mode  and  the  critical  stress  is  only 
42%  of  the  Euler  value  (lower  than  the  column  buckling  stress  for 

a  *  1500  mm).  In  the  case  of  a  panel  of  length  2670  mm,  chosen  to 
give  approximate  equality  of  local  (n  ■  7)  and  column  (n  -  1) 
buckling  stresses,  the  column  mode  is  only  slightly  Influenced  by 
local  deformation. 
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Nonlinear  and  Post-buckling  Behaviour 

14.  Nonlinear  load-deformation  relationships  have  been  computed 
for  four  panels  with  lengths  a  =■  3500,  2670,  2000  and  1000  mm, 
supplementing  results  reported  previously  (2)  for  the  experimental 
panel  L2  (a  -  3060  mm),  which  were  found  to  correlate  fairly  well 
with  test  data.  Computed  load-shortening  curves,  ie  average 

I  I 

compressive  stresses  ox  plotted  against  average  strains  ex,  are 

shown  in  Figure  5.  Mid-span  lateral  displacements  of  stiffeners 
are  shown  in  Figure  6s  positive  displacements  (at  section  AB)  are 
towards  and  negative  displacements  (at  CD)  away  from  the  stiffener 
outstand. 

15.  For  panels  of  length  3500,  2670  and  2000  mm,  compressive  loads 
reach  their  effective  maxima  at  or  just  below  initial  overall 
buckling  stress  levelss  collapse  might  be  expected  to  occur  at 
these  load  levels  as  a  result  of  outer-fibre  material  failure  at  a 
strain  of  about  0.02.  Computed  outer-fibre  stiffener  strains  ex  at 

mid-span  positions  (AB  and  CD)  are  shown  for  each  panel  in  Figure  7. 
Computed  deformation  of  the  longest  panel  at  peak  load  is  shown  in 
Figure  8. 

16.  In  the  shortest  panel  (a  =  1000  mm),  a reaches  a  plateau  at 

about  120  MPa,  substantially  higher  than  the  initial  n  =  1  buckling 
stresses  but  lower  than  the  equation  8  Euler  stress  (161  MPa):  it 
appears  that  this  load  corresponds  to  "column"  buckling  in  which 
overall  flexural  rigidity  of  the  stiffened  panel  is  reduced  by  loss 
of  effective  width  caused  by  local  buckling.  Computed  post-buckling 
deformations  of  the  panel  at  average  compressive  stresses  of  73  and 
120  MPa  are  shown  in  Figure  9.  It  is  evident  from  this  diagram  and 
from  Figure  6  that  significant  "column"  displacements  develop  at 
the  higher  load  level.  Short-wavelength  compressive  buckling  of 
the  stiffener  webs  is  also  apparent. 

17.  It  may  be  noted  from  Figure  5  that  in  the  two  shortest  panels 
(a  »  2000,  1000  mm)  occurrence  of  local  buckling  at  an  average 
stress  of  about  35  MPa  results  in  a  reduction  of  about  40%  in  the 

I  f 

compressive  stiffness  (dox/dex)  of  each  panel.  A  further  reduction 

of  about  10%  in  axial  stiffness  occurs  at  an  average  stress  of 
about  65  MPa,  corresponding  approximately  to  the  symmetric  n  »  1 
buckling  stress.  Associated  distributions  of  compressive  stress 
over  the  cross-section  and  loss  of  effective  width  are  indicated  in 
Figure  10.  If  the  flexural  rigidity  El  of  the  cross-section  is 
modified  to  reflect  this  loss  of  effective  width,  equation  8  gives 
an  Euler  buckling  stress  of  132  MPa,  about  10%  higher  than  the  load 
plateau  indicated  by  nonlinear  analysis. 

Interactive  Effects 

18.  Interactive  buckling  effects  in  a  longitudinally  stiffened  FRP 
panel  under  longitudinal  compression  may  be  summarized  as  follows. 
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(i)  Buckling  and  post-buckling  behaviour  will  usually  be 
strongly  influenced  by  coupling  between  strips  of 
laminate  forming  the  shell  and  stiffeners:  this  effect 
invalidates  the  use  of  simple  formulae  for  evaluation 
of  local  buckling  stresses  (except  to  provide  crude 
upper  and  lower  bounds)  and  can  only  be  accounted  for 
by  use  of  folded-plate ,  finite-element  or  finite-strip 
analysis.  In  long  panels  for  which  the  overall  (n  *  1) 
initial  buckling  stress  is  of  the  same  order  as  or  less 
than  the  local  buckling  stress  and  in  which  the  pre¬ 
ferred  local  buckling  wavelength  is  a  small  fraction 

of  the  span  a,  interaction  between  overall  and  local 
initial  buckling  is  slight  and  the  former  may  be 
examined  with  reasonable  accuracy  using  the  Euler 
formula.  In  short  panels,  where  the  local  initial 
buckling  stress  is  substantially  less  than  the  overall 
buckling  stress  and  the  preferred  number  of  buckling 
half-waves  (n)  is  small,  overall  instability  may  be 
strongly  influenced  by  local  deformation  and  may  occur 
at  a  stress  substantially  less  than  the  Euler  value. 

(ii)  Local  buckling  of  the  shell  laminate  between  stiffeners 
is  associated  with  loss  of  effective  width,  as  illus¬ 
trated  in  Figure  10,  causing  an  increase  in  stiffener 
stress  which  will  tend  to  accelerate  outer-fibre 
material  failure,  together  with  a  reduction  in  the 
incremental  (tangent)  flexural  rigidity  of  the  cross 
section  affecting  overall  beam-column  behaviour.  This 
interactive  effect  is  likely  to  be  most  marked  in  short 
stiffened  panels  where  local  buckling  precedes  overall 
instability.  It  is  part  of  the  reason  why  computed 
load-displacement  curves  (Figures  5  and  6)  for  panels 
with  a  =  2000  mm  and  a  =  2670  mm  reach  a  plateau,  where 
collapse  is  likely  to  occur  as  a  result  of  outer-fibre 
material  failure,  somewhat,  below  the  initial  overall 
buckling  stress.  In  the  case  of  the  shortest  panel 
examined  (a  =  1000  mm)  the  overall  (n  =  1)  buckling 
mode  is  dominated  by  local  deformation;  despite  loss  of 
effective  width  associated  with  local  buckling,  com¬ 
puted  load-displacement  curves  climb  to  a  level  well 
above  the  initial  n  -  1  buckling  stress.  It  should  be 
noted  however  (as  demonstrated  in  Reference  2)  that 
compressive  failure  of  a  short  panel  may  be  caused 
prematurely  by  debonding  of  stiffeners  induced  by  local 
buckling.  Failure  may  also  be  precipitated  by  outer- 
fibre  material  failure  of  the  locally  buckled  shell 
laminate. 

(iii)  Bending  of  stiffeners  associated  with  overall  beam- 
column  deformation  may  induce  significant  secondary 
compressive  or  tensile  stress  in  the  shell  laminate 
between  stiffeners.  Where  bending  is  towards  the 
stiffener  outstand  such  secondary  stress  is  compress¬ 
ive,  accelerating  the  occurrence  of  local  buckling, 
loss  of  effective  width  and  hence  loss  of  overall 
flexural  rigidity;  where  bending  of  stiffeners  occurs 
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away  from  the  outstand,  secondary  stress  in  the  shell 
is  tensile,  tending  to  inhibit  local  buckling  and  loss 
of  effective  width.  This  effect  is  illustrated  in 
Figure  8.  A  consequence  is  that  lateral  deformation  of 
a  multi-span  stiffened  panel  tends  to  grow  in  an 
un symmetrical  manner  with  displacements  towards  the 
stiffener  outstand  exceeding  displacements  in  the 
opposite  direction.  A  two-span  "continuous  beam" 
model,  as  shown  in  Figure  3,  is  necessary  to  account 
for  this  form  of  behaviour  in  a  numerical  analysis. 

The  effect  is  however  found  to  be  small  for  cases 
examined  in  the  present  papers  over  the  range  of 
deformation  shown  in  Figure  6  the  maximum  difference 
between  "upward"  and  "downward"  displacements  is  10%. 

(iv)  Approximate  equality  of  critical  stresses  referring  to 
local  (n  =  7)  and  column  (n  =  1)  buckling  modes  (for 
the  panel  of  length  2670  mm)  does  not  appear  to  cause 
any  special  enhancement  of  the  interaction  effects 
discussed  above. 

DESIGN  IMPLICATIONS 


19.  Current  design  recommendations  for  GRP  ships  (2)  specify  that 
factors  of  safety  of  2.0  should  be  maintained  against  local  and 
overall  buckling  of  longitudinally  stiffened  panels  under  compress¬ 
ive  load,  reduced  to  1.5  in  the  case  of  local  buckling  where  posi¬ 
tive  measures,  eg  use  of  bolts  or  resilient  adhesive  with  high  peel 
strength,  are  employed  to  enhance  stiffener  attachment.  As  dis¬ 
cussed  in  Reference  2,  these  factors  are  intended  to  account  for 
typical  initial  imperfections  but  not  for  creep  or  for  variability 
and  long-term  degradation  of  material  properties,  which  should  be 
alloved  for  separately  by  application  of  supplementary  partial 
safety  factors. 

20.  Various  forms  of  interaction  between  local  and  overall  buck¬ 
ling  have  been  identified  in  the  present  numerical  study.  For 
panels  with  cross-sectional  geometry  similar  to  that  of  Figure  2  it 
appears  that  interaction  effects  are  accounted  for  adequately  by 
the  margins  indicated  above  and  that  no  additional  partial  safety 
factor  is  necessary,  even  for  the  case  where  critical  stresses 
referring  to  local  and  overall  buckling  modes  are  approximately 
equal.  It  seems  likely  that  these  conclusions  will  apply  equally 
to  most  practical  stiffened  panels:  scope  clearly  exists,  however, 
for  further  exploration  of  interactive  buckling  effects  in  panels 
of  differing  geometry  and  for  different  patterns  of  initial 
deformation. 
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FIG  1  Longitudinally  Compressed  GRP  Panel 
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FIG  2  Panel  Geometry  and  Material  Properties 
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Selected  Initial  Buckling  Modes  Computed  by  Folded  Plate  Analysis 
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FIG  5  load-Shorteninij  turves  tor  Stitferej  » anels 


FIG  6 


Mid-Span  Lateral  Displacement  of  Stiffeners 
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FIG  10  Compressive  Stress  Distributions  Across  Shell  Laminate  Showing 
Reduced  Effective  Width 
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